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Abstract 

Ultra-thin magnetic tunnel junctions with low resistive MgO tunnel barriers are prepared to ex- 

■^ 

t— 5 amine their stability under large current stress. The devices show magnetoresistance ratios of up to 

00 

psj 110 % and an area resistance product of down to 4.4 rj/xm . If a large current is applied, a reversible 

f"! resistance change is observed, which can be attributed to two different processes during stressing 

"o 

^ and one relaxation process afterwards. Here, we analyze the time dependence of the resistance 

^ and use a simple model to explain the observed behavior. The explanation is furter supported by 

g 
. numerical fits to the data in order to quantify the timescales of the involved phenomena. 
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I. INTRODUCTION 

The switching of a ferromagnet's magnetization by injection of a spin polarized current 
(Spin Transfer Torque STT switching) is an important new phenomenon in magnetism PQ [2] . 
It enables applications such as spin torque nano-oscillators (STNO), Magnetic Random 
Access Memory (MRAM), programmable magnetic logic and sensors [3] using the tunneling 
magnetoresistance (TMR) effect [H El [6] . The conventional switching scheme of field pulses 
generated by current lines requires large currents and is not scalable, which led to alternatives 
such as heat assisted [7j or the already mentioned STT switching [S]. In state-of-the-art 
MgO based magnetic tunnel junctions (MTJs), STT requires current densities of about 
1 X lO^A/cm^ [9]. This leads to a bias voltage of about 100 mV, if the area- resistance 
products RA are in the order of lOQfim^. 

Individual STNOs are strongly limited in their output power. This deficiency has to 
be overcome in order to exploit the benefits of current-tunable narrow-linewidth microwave 
generation for future application. Although an array of coupled oscillators may offer a 
path to higher power output (TU], the individual STNOs are envisioned to be MTJ-based 
and operating close to the degradation current limit [Sj [H]. Please note that STNOs are 
operated at direct currents in contrast to pulsed currents in the case of MRAMs. 

A constant trend to reach lower area resistance products forced the fabrication of MTJs 
with barriers as thin as 10 A. Further reduction of the tunnel barrier thickness is becoming 
more and more challenging. Recently, reversible resistance changes and an atypical tem- 
perature dependence of the resistance were observed in thin AlOx MTJs when reaching 
current densities of 10^ - lO^A/cm^ [12]. These resistance changes do not depend on the 
relative orientation of the magnetic layers [12] and are attributed to electromigration in 
nano-constrictions of the insulating barrier [13l[Tlj. Two opposite relaxation processes were 
found suggesting two independent mechanisms acting simultaneously inside the MTJs [TB] . 

Here, we study the effect of current stressing on ultra-thin MgO barriers. In these devices, 
resistance changes take place at a current density of about 10^ A/cm^. This shows that fur- 
ther reduction of the STT switching currents is crucial in order to achieve reliable operation 
of STT-MTJ based devices. We distinguish two effects during the current stress and one 
relaxation process. Furthermore, numerical fits to the experimental data are presented. 



II. RESULTS 

The TMR systems with low resistive MgO barriers is sputter deposited in a Singulus 
TIMARIS II tool. The film sequence is Ta 3/Cu-N 90/Ta 5/Pt-Mn 20/Co-Fe 2/Ru 0.75/ 
Co-Fe-B 2/MgO 1.3/Co-Fe-B 3/Ta 10/Cu 30/Ru 7 (all values in nm). The composition of 
the compounds is: PtsyMugs, CoyoFeso, and Co66Fe22Bi2. Elliptical TMR elements with sizes 
in the range of 0.018 yum^ to 0.095 /im^ are prepared by e-beam lithography in combination 
with ion beam etching. By applying different Ar pressures during MgO deposition two 
sample sets with an area-resistance product of AAQfim^ (set 1) and 9.7Qfim^ (set 2) are 
produced. 

All measurements are performed using a constant voltage source at room temperature. 
By applying a voltage in the range of 200 mV to 800 mV the samples are stressed by a 
current density of 1 - lOx 10^ A/cm^. A magnetic field of ±400 Oe provides a stable parallel/ 
anti-parallel orientation of the magnetic layers during the measurements. The resistance is 
monitored at a read-out voltage of 20 mV with a rate of 1 Hz in order to overcome non- 
linear J(y)-contributions and thermal drifts related with heating processes (heating rates 
due to current load have typical time constants well below Ins [15]). Each measurement is 
performed on an individual element which was not stressed before. This is necessary since 
the properties of the barrier may change irreversibly by the stressing. 

Figure [I] shows the MTJ resistance upon stressing with positive voltage polarity (electrons 
flow from the bottom to the top electrode). The area-resistance product is RA = 9.7r2yum^. 
The inset shows the magnetic minor loops before (black) and after (gray) stressing exhibiting 
TMR ratios of the junction of 104 % and 108 % at 20 mV, respectively. This increase in TMR 
gives evidence, that the barrier itself does not degrade during current stressing. As the device 
is stressed by j = +1.5 x 10^ A/cm^ one can see a decay of the resistance with time, although 
the films had been annealed for 1.5 h at 360°C prior to patterning. The initial resistance 
of -Rini = 465 Q is reduced reaching a saturation resistance of Rgat = 442 Q within 30 min. 
Numerically, the measured time dependence of the resistance can be fitted with moderate 
quality by an exponential function: 

R{t) = Rs.t + Ai? ■ exp i-t/r) (1) 

with AR = Rini — Rsat > 0. This would represent a process characterized by a decay time 
of r = 3.15 min and leads to the dashed line in Fig. [T]as will be further discussed below. 



Stressing the MTJs with currents of the negative polarity leads to a substantially different 
response shown in Fig. |2} First, the resistance rapidly increases from Ri^^ = 82.5 i7 to 
-Rmax = 84.5 Q, followed by a slow decay. No saturation of the resistance can be seen within 
the measured time range. 

Thus at least two processes are involved: A fast, polarity dependent resistance change 
followed by a slow resistance decrease for both polarities: 

R{t) = Rs.t + Ai?i exp(-t/ri) + Ai?2 exp(-t/r2) (2) 

with ARi > (< 0) for positive (negative) polarity and AR2 > for both polarities. These 
two exponential decays of different sign give a very good fit to the data of Fig. [2] (solid line) 
with Ti = 2.8 min ± 0.2 min and T2 = 15.1 min ± 2.7min. Looking back to Fig. [T]the fast 
process might be hidden in the decrease because the sign depends on the polarity of the 
applied voltage: Two exponential decreases would lie on top of each other. In fact with Eq. 
[2] a better fit to the data of Fig. II] can be obtained (solid line). In this case, however, the fit 
parameters ti and T2 are strongly correlated and can not be used for a quantitative analysis. 

The dependence of the time constant r (Fig. [I| dashed line) on the RA product of the 
junctions is shown in Fig. |3]for different positive stressing voltages. For increasing voltage a 
decrease in r is observed. At the same time, for a given voltage, r increases with increasing 
RA product. 

Figure |4] shows the current induced resistance variation AR as a function of the current 
density. A decrease of the resistance can be seen reaching 13% at a current density of 
j = 4.5 X lO^A/cm^ for sample set 2. For sample set 1 AR is reduced by a factor of 
approximately 0.45 due to the lower RA. This can be seen from Fig. |4] where the black line 
shows a linear fit to the data of sample 2 and the green line is generated by changing the 
slope of the black line by the RA ratio which is 8.9Qfim^ to 19.8Qfim^ for the antiparallel 
state. In contrast to this, the magnetic orientation of the electrodes (compare green and 
orange squares) produces no significant differences. This leads to the conclusion that the 
current density and not the electric field is the main parameter that determines AR. 

If we investigate the relaxation of the resistance, i. e. the resistance change after the 
current stress is released, one can reveal a third process involved. The red curve in Fig. [5] 
is equivalent to the plot shown in Fig. [l| the blue curve is equivalent to the plot in Fig. |2} 
The green curve in Fig. [s] illustrates the relaxation process: If no large voltage/ current is 



applied (20niV), the resistance slowly reaches the initial level i?o- This can also be fitted 
with an exponential function with typical relaxation times of 100 min. 

III. DISCUSSION 

For the interpretation of the data, one has to consider that the resistance of a (magnetic) 
tunnel junction is governed by the barrier thickness d and the barrier height $. This re- 
lationship can be described by Brinkman's equation [16]. For barrier height (0.5 - 3.5 eV) 
and thickness (1 - 1.5 nm) that correspond to our junctions, the Brinkman equation can be 
reduced to 

i?(ci,$) = ci-exp(c2rf$^/^), (3) 

where Ci and C2 are constants. Thus, the resistance changes as given by Eq. [T] could be 
explained by a variation of d as well as $. 

If we assume that the effective barrier thickness changes by a certain fraction with a nom- 
inal thickness of Ac/, Brinkmans's equation leads to Ad ~ 0.1 A, if the resistance change 
is about 5%. For a change of $, a resistance change of 5 % leads to A$ k, lOmeV. Thus 
very small changes in the barrier characteristics can lead to the observed resistance changes. 
Provided small changes of c? or $ results in a linear approximation of Eq. |3] and the exper- 
imental findings (Eq. [I]) lead to an exponential function for either d{t) oc (1 + exp(— t/r)) 
or $(t) oc (1 + exp(— t/r)) as well. 

Since the layer stack is symmetric (CoFeB/ MgO/ CoFeB), the same resistance change 
for positive and negative bias is expected. The symmetry break is induced by the sequence 
of preparation. The lower ferromagnet is more likely to be contaminated with oxygen |17] , 
because the MgO is deposited on top of the ferromagnet by rf-sputtering from an MgO 
target. This oxygen on top of the ferromagnet can act as an additional barrier [TH] and also 
lead to oxygen deficiencies in the lower part of the tunnel barrier. Also the effective barrier 
height $ is modified by the presence of oxygen vacancies [T9| [20] . 

Therefore, a possible explanation for the observed behavior is a small displacement of 
these oxygen vacancies/ defect atoms driven (fast) in the direction of the electron flow 
increasing or decreasing the effective thickness of this layer (electromigration) [14J and (slow) 
by thermal activation. These processes are reversible since the atoms can slowly fall back 
into their initial positions of lower energy once the current is switched off. 



IV. SUMMARY 

In summary, we investigate magnetic tunnel junctions with MgO barriers with magneto- 
resistance ratios of up to 110 % and area-resistances of around lOQfim^. The low resistance 
of these devices allow us to drive large currents trough the structures, which is necessary 
for use in spin transfer torque memory or logic cells. We observe three reversible resistance 
changes of the junctions: One fast, polarity dependent term that could be associated with 
electromigration. One slower, polarity independent term active during current stress, that 
is suggested to be of thermal origin. Finally, a very slow relaxation takes place when the 
current stress is released, which also can be explained by thermally activated diffusion. 
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FIG. 1: Resistance as a function of time for the positive current polarity. A fast resistance decrease 
can be observed during the first 5 min. The data is fitted using Eq. [T] and Eq. |2] as shown by the 
dashed and sohd hne, respectively. The inset shows the minor loop of the investigated MTJ before 
(black) and after (gray) stressing. 
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FIG. 2: Resistance as a function of time for the negative polarity. A fast increase of the resistance 
is followed by a slow decrease. The data is fitted by Eq. [2] as shown by the solid line. 




FIG. 3: Decay time r presented by a color scale as a function of voltage and area-resistance. A 
decrease with increasing voltage is observed. Increasing the area-resistance at a given voltage leads 
to a larger r. 
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■ set 1 AP (8.9ilnm^} • 
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FIG. 4: Resistance change as a function of current density. The black line represents a linear 
fit to the results obtained from sample set 2. The slope of the green line is calculated from the 
difference in RA between set 1 and 2. Parallel and antiparallel alignment of the magnetic electrodes 
is indicated by P and AP, respectively. 
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FIG. 5: Resistance behavior under stressing with positive polarity (red), during the relaxation 
(green), stressing with negative polarity (blue), and relaxation (green). The resistance Rmi as 
obtained prior to the stressing is indicated by the dashed line. 
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